Physical Properties of the B and Be Star 
Populations of h and \ Persei 



Amber N. Marsh Boyer 1 , M. Virginia McSwain 1,2 , Christina Aragona 1 , 

and Benjamin Ou-Yang 3 

Department of Physics, Lehigh University, 16 Memorial Drive East, Bethlehem, PA 18015 
anm506@lehigh.edu, mcswain@lehigh.edu, cha206@lehigh.edu, ouyang@chara.gsu.edu 



Of 



ABSTRACT 

We present a study of the B and Be star populations of the Double Cluster h and \ Persei. Blue 
optical spectroscopy is used to first measure projected rotational velocity, V sin i, effective surface 
temperature, T c ff, and surface gravity, log g, for B-type sample stars, while available Stromgrcn 
photometry is used to calculate T e s and log g for the Be stars showing emission. In our sample of 
104 objects for which we measured these stellar parameters, 28 are known or proposed Be stars. Of 
these Be stars, 22 show evidence of emission at the times of our observations, and furthermore, we 
find evidence in our data and the literature for at least 8 transient Be stars in the clusters. We find 
that the Be stars are not rotating near their critical velocity, contrary to the results of studies of 
similar open clusters. We compare the results of our analysis with other previous studies, and find 
that the cluster members are more evolved than found by Huang & Gies but still retain much of 
their initial rotational angular momentum. 
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NGC 869 and NGC 884 (h and X Persei, respec- 
tively) are a well known double open cluster, visi- 
ble in the Northern hemisphere, and have been the 
focus of many studies over the years. The early 
1900's saw a number of studies attempting to deter- 
mine cl uster membership, positions, and radial ve- 
locities dMessowlll913t lAdams fc VanMaanenlll913t 
Hertzsprunj 1922 ). By the 1960' s more extensive 
studies, such as that of ISlettebakl (Il968l) . were be- 
ing conducted to determine spectral types for the 
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cluster constituents. More re cently an ex t ensive 
study has been conducted by ICurrie et all (|2f)10h 
in which they investigated the general properties 
and membership of the cl usters. Their res u lts, in 
agr eement with those o f iBragg fc Kenvonl (|2005 ) 
and lSlesnick et aL ( 2002 ). find that the clusters are 
incredibly similar, having common ages of ~13-14 
Myrs, distance moduli dM = 11.8 - 11.85 (~ 2, 200 
pc), and reddenings of E(B-V) ~ 0.52-0.55. They 
also estimate a total mass of at least 20,000 M Q for 
the clusters. 

One of the prominent motivations for our study is 
that these young open clusters are rich in Be stars. 
As early a s the 1920 ' s, obs ervational studies con- 
ducted by iTrumplerl (|l926l) and others noted the 
presence of emission in the hydrogen lines of many 
of the brightest B-type cluster members. Modern 
studies of the cluster have shown that upwards of 
30% of t he brightest B-typ e stars are known to be 
Be stars ( Keller et al.l 2001 ) . In a study conducted 
with Spitzer, ICurrie et al. f (l2008h investigated the 
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lower mass stellar population for mid-infrared ex- 
cesses due to the presence of protoplanetary disks. 
They also identified 57 Be stars and candidates ex- 
hibiting excess emission at 24/im, which helped to 
motivate our study to follow up these candidates 
and confirm their Be nature. Of their stars, 21 had 
previously been identified as showing emission, and 
20 of their stars are included in the present study. 

The modern working definition of a Be star is 
given as "a non-supergiant B star whose spectrum 
has or had at some time, one or mo re Balmer lines in 
emission" ([Porter fe Riviniusll2003l ) . While this gen- 
eral definition also encompasses objects such as the 
well known Algol binary systems and Herbig Ac/Be 
stars, classical Be stars are further delineated as 
having circumstellar line emission formed in an opti- 
cally thin equatorial disk, l ow-order line profile var i- 
ations, and rapid rotation (jPorter fe Riviniusll2003l) . 
These Be disks are comprised of warm gaseous mate- 
rial ejected from the stellar surface during outburst 
events. The gas is then pulled into a gravitationally 
bound orbit about the stellar equator. 

It is well established that as a population, Be 
stars rotate faster than than their non-emission, B- 
type counterparts. Precisely why this is the case, 
however, is still debated. There are three primary 
theories as to why Be stars are rapid rotators: they 
may have been born as rapid rotators, spun up by 
mass transfer in a close binary system, or spun up 
during the main sequence evolution of B-type stars. 
The observed rotation rat es of Be stars are > 60- 80% 
of their critical velocity ( McSwain et al. 20081) . at 
which point the gravitational and centrifugal forces 
are balanced, although recent results suggest that 
this t hreshold may be mass dependent (|Huang et al.l 
2010h . The main sequence lifetimes of these objects 
are likely extended as a direct result of their rapid 
rotation, as this fosters rotational mixing of their 
stellar interiors and reple nishes their hydrogen cores 
(jMevnet fe Maederllioool ). However, rapid rotation 
alone is not enough to spur the photospheric mate- 
rial of these stars to form the disk structures they 
host. It is likely that other weaker processes, such as 
non-radial pulsations (NRPs), are needed to provide 
the additional angular momentum ne cessary for this 



material to leave the stellar surface ( Riviniu s et al 
20011: IPorter fe Riviniusll2003t iMcSwain et al.ll2008t 



B-type and Be stars have become important tar- 
gets for asteroseismology. In hot stars, different pul- 
sation frequencies provide information about differ- 
ent layers of the stellar interior. The information 
gleaned from studies targeting B stars with NRPs 
is being used to improve current stellar structure 
and stellar evolution models. Doing so, however, re- 
quires accurately determined stellar surface parame- 
ters in order to set appropriate boundary conditions 
for these models. Accurate measures of effective 
temperature and surface gravity are particularly es- 
sential to determining stellar radii for these models, 
as well as for determining stellar ages and evolu- 
tionary spin-down. Many B-type and Be stars in h 
and x P er have been found to host NRPs, and cur- 
rently there are on-going campaigns to observe and 
characterize the natur e of the variable pulsations 
found in these stars ferzesiriski fc Pigulskil Il997t 
Krzesihski et al.l 119991 Isaesen et al.ll2010h . The in- 



terest in the pulsating stars of h and x P er necessi- 
tates improved measurements of stellar parameters 
for the cluster. 

There is an on-going debate in the massive star 
community regarding the evolution of angular mo- 
mentum of B-type stars. With their abundance of 
B-type stars, h and x P er are two of t he many stel 



lar clusters at the center of this debate. Strom et al 



(2005) find that the present-day rotation rates of 
these stars are set by environmental characteristics 
of the natal clouds in which they formed, with little 
change o ver the main se q uence s tellar lifetime. The 
work of [Huang fc Gied (<2006ah and iHuang et al 
(2010), lowever, indicates that the observed rota- 



Cranmerl l2009) . A growing number o f Be stars have 



been identified to exhibi t NRPs (see iRivinius et al 
2003 ; ; lEmilio et alll2010l ) . 



tion rates of B-type stars are due less to the initial 
birth-line rotational rates of the stars and more to 
evolutionary spin-down or mass transfer in binaries. 
Both studies also observe that B-type stars in clus- 
ters have, on average, significantly higher rotation 
rates than field B-type stars. Resolving the angu- 
lar momentum problem also requires a new exam- 
ination of the stellar parameters and evolutionary 
states. 

In this work we present stellar properties de- 
termined from blue optical spectroscopy of a sam- 
ple of B and Be stars in the clusters. Section 2 
provides details of spectroscopic observations taken 
with various instruments. In Section 3 we discuss 
our methods for determining V sin i, T e g, and log g 
for both the B and Be star populations. We feel 
that in order to study Be stars relative to their 
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non-emission peers it is important to apply a uni- 
form technique for B and Be star analysis (in as 
much as possible). With this homogenous and in- 
dependently determined dataset we will be able to 
further investigate the nature of the clusters' Be 
star population. Section 4 quantitative l y com pares 



our results w i th tho se of Huang fc Gied 



_ 2006aj) and 
Strom et al. (2005), who use different methods to 
determine the same physical parameters. We dis- 
cuss the impact of our disagreements in the context 
of the evolution of stellar angular momentum. Fi- 
nally, in Section 5 we draw our conclusions regarding 
the clusters' massive star populations. We plan to 
use the results of this work to examine the spectral 
energy distributions (SEDs) of the star + disk sys- 
tems among the Be star population in a forthcoming 
paper. 

2. Observations 

We have obtained spectra for a total of 104 mem- 
bers of NGC 869 and NGC 884 during multiple ob- 
serving runs: 2005 November using the Kitt Peak 
National Observatory (KPNO) Wisconsin Indiana 
Yale NOAO (WIYN) 3.5 m telescope with the Hy- 
dra multifiber spectrograph; 2010 August using the 
Wyoming Infrared Observatory (WIRO) 2.3 m tele- 
scope with the Long Slit spectrograph; 2011 Novem- 
ber using the KPNO 2.1 m telescope with the Gold- 
Cam spectrograph; and 2012 January using the 
0.9 m KPNO Coude Feed (CF) telescope with the 
Coude spectrograph. The UT dates, wavelength 
range, resolving power, number of targets, and in- 
strumental setup details for all runs are summarized 
in Table [U 

All of the spectra obtained at the WIYN 3.5 m 
with the Hydra spectrograph have been zero cor- 
rected using standard routines in IRAP0, and have 
been flat-fielded, wavelength-calibrated, and sky- 
subtracted in IRAF using the dohydra routine. 
The Hydra observations obtained by M. Virginia 
McSwain in 2005 consist of 7 exposures of NGC 869, 
totaling 2.25 hrs, and 5 exposures of NGC 884, to- 
taling 2 hrs. Each exposure has been wavelength 
calibrated via a CuAr comparison spectrum both 
before and after the cluster observations. For each 
of the two configurations, the exposures have been 
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transformed to a common heliocentric wavelength 
grid and co-added to produce good signal-to-noise 
for each star. The spectra were rectified to a unit 
continuum by fitting line-free regions. 

Before any reduction or calibration routines were 
applied to the WIRO 2.3 m data, all spectra (ob- 
ject and comparison) were corrected for bit-flip er- 
rors with the rfits routine in IRAF. A CuAr cali- 
bration lamp source was used to obtain wavelength 
calibration spectra before and after every object 
spectrum. The spectra were then zero-corrected, 
flat-fielded, wavelength-calibrated, and rectified to 
a unit continuum using standard slit spectra rou- 
tines in IRAF. Fainter objects that required mul- 
tiple exposures were co-added prior to continuum 
rectification to improve signal-to-noise. Given the 
poorer resolution of blue WIRO spectra, we cannot 
use these data with our spectral fitting techniques 
described in Section 3. We therefore only use the 
blue and red WIRO spectra to identify Be stars ex- 
hibiting emission at the time of our observations. 

Spectra from the KPNO 2.1 m instrument 
have been similarly zero-corrected, flat-fielded, and 
wavelength-calibrated using the standard routines 
found in IRAF. Once wavelength calibrated via the 
HeNeAr comparison lamp spectra, which were taken 
before and after every object spectrum, the data 
were rectified to a unit continuum. 

In a similar manner, the spectra obtained with 
the 0.9 m KPNO CF telescope have been zero- 
corrected, flat-fielded, and wavelength-calibrated 
using the standard routines in IRAF. ThAr com- 
parison spectra were taken every one to two hours 
during the run. The data were then rectified to a 
unit continuum. 

3. Physical Parameters from Spectral Mod- 
els 

The first component of our population study 
is to determine basic parameters for the clus- 
ter constituent s using a methodology devised by 
McSwain et al. (120081) . Using ground based opti- 
cal spectroscopy (4000-5200 A) and model fitting 
techniques, we can determine V sin i, T c ff, and log 
g for each star. To obtain these measurements, we 
compare our observed spectra to grids of model B- 
type stars, determining a best fit to the data by 
minimizing the mean square of the deviations rms 2 . 

We begin our analysis by making a rough esti- 
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mate of T e g- and log g for a star, and then we com- 
pare the He I AA4387, 4471, 4713, and Mg IIA4481 
lines with the Kurucz ATLAS9 models ( Kurucd 
19941 ) to determine V sin i. These lines are used 



can determine the critical velocity 



because their broadening is dominated by the ef- 
fects of rotation, thus yielding a better indication of 
V sin i. We then take a weighted average of the four 
values determined from each of the lines to give our 
measured value of V sin i. The error, Ay sin i, is 
determined by the offset from the measured best-fit 
value that increases the rms 2 by 2.7rms 2 /iV. Here, 
N is the number of wavelength points within the fit 
region. A sample fit determined for He I A4387 in 
NGC 869-90 is shown in the bottom panel of Figure 
[TJ Our results for V sin i and its errors are listed in 
columns 2 and 3 of Tables [2] and [3] 

In Be stars it is possible that the He I lines may 
contain weak emission from the circumstellar disk, 
which may partially fill the absorption features and 
narrow the overall line profile. So while we have 
made measurements of V sin i where we can for Be 
stars in our sample, these values should be consid- 
ered as lower limits. 

We measure a mean V sin i =157 km s _1 with 
a standard deviation of 89 km s _1 for the normal 
B-type stars of both clusters, including binary sys- 
tems. Assuming an average inclination angle of 
i = 60°, this gives a mean V cq — 181 km s _1 for 
our sample of B stars. For the Be stars we measure 
a mean V sin i — 205 ± 81 km s _1 , with a mean 
V cq = 237 km s _1 . From this it is clear that the 
Be stars in these clusters are, on average, rotating 
somewhat more rapidly than their B-type counter- 
parts. 

The cumulative V sin i distributions for all Be 
and normal B-type stars in both NGC 869 and 
NGC 884 are shown in Figur e [H In comparison to 
other young open clusters (see McSwain et al. 20081 
2009), we find that the Be star population of NGC 



869 and NGC 884 are rotating surprisingly more 
slowly that expected in comparison to their B-type 
counterparts. Using the two-sided Kolmogorov- 
Smirnov (K-S) statistical test, we investigate the 
null hypothesis that the distributions of B-type and 
Be stars differ. The K-S test indicates a 7.6% chance 
that the two populations are drawn from the same 
sample. Using the values of mass, M*, and radius, 
R+, discussed below for all stars in our sample, we 



Vc.rit — 



Re 



(1) 



for our stars. For simplicity in this expression, we 
assume that the polar radius of the star, R p , is equal 
to R+, and that a rotationally distorted star has an 
equatorial radius R c — 1.5i? p . With this, we find 
a mean V CT it of 430 km s _1 for the B-type and Be 
stars in these clusters. 

Having determined values of V sin i for each star, 
we turn again to model spectral fitting to determine 
values for T e ff and log g from the H7 line at 4340A. 
The hydrogen Balmer lines are particularly sensitive 
to T c ft and log g, making them ideal for determining 
these quantities accurately. The method outlined 
here was employed for normal B-type stars and Be 
stars with no emission present in this work. For 
stars having T e s < 15000 K we employ the meth- 
ods of iHuang fe Giesl (|2006ah . who use H7 line pro- 
files generated by the line-blanketed, local thermo- 
dynamic equilibrium (LTE) Kurucz ATLAS9 and 
SYNSPEC codes. The "virtual star" models pro- 
duced by their code simulate spherically symmet- 
ric stars with constant T e g and log g across their 
surface, and the spectra of these model stars are 
then used to measure these values and their er- 
rors of our observed spectrum, similar to our pro- 
cedure for V sin i. The errors AT e g and Alog g 
are determined from the quadratic sum of the V 
sin i propagated errors and the errors due to the 
intrinsic noise in the observed spectrum. For hot- 
ter stars, LTE models should systematically under- 
estimate T e ff as non-LTE effects alter the equiva- 
lent width of the H7 line we are measuring. Hence 
for stars having T s > 15000 K, we use instead 
the metal line-blanketed, non-LTE, plane-parallel, 
hydrostatic TLUSTY BSTAR2006 model spectra 
(jLanz fc Hubenvll2007h . A sample fit of the H7 line 
in NGC 869-90 is shown in the top panel of Fig- 
ured The errors AT ff and Alog g are determined 
from the values that produce an rms 2 no more than 
2.7rms 2 /iV greater than the minimum rms 2 . Our 
results for T e g, log g, and their respective errors are 
shown in columns 4-7 of Table [2] 

29% of the stars in our sample arc rapid rota- 
tors, having measured values of V sin i in excess 
of 200 km s _1 . At such significant rotational ve- 
locities, assuming a spherical shape for these stars 
is no longer plausible given the substantial centrifu- 
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gal forces distorting the stars into oblate spheroids. 
This rotational distortion produces significant dif- 
ferences in both the temperatures and surface grav- 
ities at the polar and equatorial regions. As the 
measured values of T e g and log g are averages across 
the visible stellar hemisphere, these rotational ef- 
fects produce lower values than expected, a phe- 
nomenon known as gravitational darkening. As the 
rotation rate for a star approaches its critical veloc- 
ity, the equatorial radius may increase by as much 
as 50%, while the polar radius remains unchanged. 
For these reasons, we conv ert our measured log q 
to log 5 P oiar, as detailed in Ifiuang fc Gied ( 2006ah . 
The authors produce detailed spectroscopic models 
to investigate the effects of such rotational distor- 
tions and determine a statistical correction factor 
for log g. This factor is averaged over all possible 
values of inclination angle, i, for a variety of stellar 
models. By bilinearly interpolating between their 
models, we converted our measured value of log g 
to log (? P oiar- This value of log g po iar is a better mea- 
sure of the true surface gravity of the star, untainted 
by the effects of rapid rotation, and provides us a 
more accurate means of comparison between slowly 
and rapidly rotating stars. The log <? po iar correc- 
tion is therefore most significant for stars rotating 
more rapidly than 50% of V CI ft. We assume that 
this conversion produces a negligible change to the 
value determined previously for Alog g. Our derived 
values of log g po iar are listed in column 8 of Tabled 
Once we have measured parameters for all of 
the B-type stars, we can then determine M* and 
_R* for each of them by interpolating values from 



Schaller et al 



the 
ary tracks 



([19921 ) non-rotating evolution- 
consistent with the slow rotation of 
most of our targets. These evolutionary tracks are 
shown plotted with T e g and log g po Ur in Figure 
[5j The errors AM* and Ai?* correspond to our 
measured AT ff and Alog g. Additionally, we have 
compared our r esults with the rotating models of 



Ekstrom et al.l (|2012h . and we find agreement be- 
tween the models to within 10%. The resulting 
values of M*, i?*, and their respective errors are 
listed in columns 9-12 of Tables [5] and |3] We have 
also checked the accuracy of our results by com- 
paring the TLUSTY BSTAR2006 model SEDs with 
the observed SEDs for our B-type sample stars, and 
find excellent agreement between our derivations of 
distance and redde ning with the accepted values of 



Currie et al. (20Jj0j). These results will be discussed 



further in a forthcoming paper. 

We note that star NGC 869-566 did not show any 
signs of emission in our initial observations, but has 
since exhibited progressively stronger emission in 
our more recent 2010 and 2012 observations; hence 
we include it among the non-emission Be stars and 
have measured T c g and log g from the H7 line it 
its 2005 blue spectrum. For the other Be stars in 
our sample we cannot determine accurate values for 
T e ff and log g via the same model fitting technique, 
as hydrogen absorption line profiles arc altered by 
emission during disk phases of these stars. Instead 
we can use Stromgren photometry available from 
the WEBDA0 database to correlate T e ff and log g for 
all of the B-type and Be stars in our sample. Non- 
emission B stars with both available Stromgren pho- 
tometry and spectral model fits were included as cal- 
ibration stars for the photometric technique. To this 



we ad d eight additional stars from lNapiwotzki et al 
(119931 ) with well-known T c g and available Stromgren 



data. We use ubvy magnitudes to first determine the 
Stromgren indices mi and c±, given by 



mi = (v -b) - (b-y) 
ci = (u — v) — (v — b). 



(2) 
(3) 



The calculation of T e ff determined by lBaloiia ( 19841 ) 
requires us to convert the c\ index to the dereddened 
index c via the expressions 

E{b-y)= E(B — V)x 0.754 (4) 
c = ci- Q.19[E(b -y)} (5) 
c= log(c +0.200), (6) 

which use the reddening values to the clusters, 
E(B - V) = 0.55 and E(B - V) = 0.52 for 



h and \ Per, respe ctively (jCurrie et al.l l2010t 
iBragg fc Kenvonll2005l) . Using these indices and 



the H/3 line magnitude (/3), we can then calculate 
T e g via the relationship 



logTeaiona = 3.9036 - 0.4816(c) - 0.5290(/3) 
- 0.1260(c) 2 + 0.0924(/3)(c) - 0.4013(/3) 2 

given by iBalonal (|l984T ) . 



(7) 



McSwain et al.l (|2008) found that this calculated 



value, TBaiona, slightly underestimates the true T e g 
of the B-type stars. Thus, we performed a linear 



5 Available online at www.univie.ac.at/webda and maintained 
by Ernst Paunzcn 
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fit to the data, shown in Figure [3j and determined 
a correction factor that will bring the two indepen- 
dent measurements into agreement. In this way, we 
are able to use the B-type stars measured by both 
methods as a calibration to yield values of T e s for 
the otherwise immeasurable Be stars. These T e ff 
and our calculated errors are listed in columns 4 
and 5 of Table El 



To determine log g, iBalonal (|1984[) advocate us- 
ing /3 and cq. Given that /3 serves as an indicator 
of log g via spectra line width, and cq is an indica- 
tor of temperature in hot stars, the relationship be- 
tween these two values can serve as a Hertzsprung- 
Russell diagram for the cluster providing a means 
for investigating the temperature and evolutionary 
trends of the stellar population. In more evolved 
giant or supergiant stars (luminosity classes III and 
I, respectively), the hydrogen lines are very narrow 
due to the lower densities in the outer atmospheres 
of these stars, decreasing the collisional rates that 
produce the pressure broadening mechanism. Main 
sequence stars (luminosity class V), which host more 
dense atmospheres and thus higher collisional rates 
and higher pressures, have broader hydrogen lines. 
In general, we do see that the values of Co and 
(3 for the B-type calibrators shown in Figure 0] 
ag ree with the relations for c lass V and III stars 
of Balona fc Shobbrookl ( 19841 ). However, the cir- 
cumstcllar disks present in Be stars (shown as filled 
diamonds in Figure |4j will artificially brighten the /3 
magnitudes of these stars, contaminating the Co — /3 
relation for these stars and our calculated values of 
log g. The spread in log g of the B-type calibration 
stars prevents us from simply applying either of the 
Co — /? relations shown in Figure |4] to the population. 
Instead, we perform a linear fit to the calibration 
star data and obtain a corrected value of (3 which 
we then use to correct our calculated Teaiona and 
finalize our calibrated, calculated temperatures and 
surface gravities for the Be stars exhibiting emission 
in our data. Since the clusters are approximately the 
same age, this single fit is appropriate. Additional 
details regarding the method to determine T e g, log 
g, M*, R+, and their respective error s for t he Be 

~ (|2008h . 



stars can be found in McSwain et al 



We 

do not perform any further correction to obtain log 
5poiar for the Be stars measured with this technique, 
given the large scatter between the calibrators' log 
Spoiar and their Strdmgren log g. The final results 
for log g and Alog g of the Be stars are listed in 



columns 6-7 of Table El 

We highlight our results regarding the Be star 
population in Table El and we note that these Be 
stars were selected from the literature based on pre- 
vious observations of spectral line emission or via 
photometric surveys needing spectroscopic confir- 
mation. Presented within Table El are a total of 
28 known or proposed Be stars within our sam- 
ple of 104 stars. For these objects we were able 
to examine the state of emission from our obser- 
vations, and then determine their stellar parame- 
ters by our spectral modeling or Stromgren pho- 
tometry methods as detailed earlier in this sec- 
tion. The broad-band Stromgren indices used to 
derive their T c g and log g are not likely to be 
affected by their rapid rotation. We do see evi- 
dence of emission in 22 of these Be stars. Inter- 
estingly, we see evide nce in our data for at least 8 
"transient" Be stars (|McSwain et all 120081 ). The 
stars NGC 869-146, NGC 869-717, NGC 869- 
1268, and N GC 884-2262 were observed as Be stars 
in the past (ISlettebakl Il985t iFabregat et al.| [l994; 



Keller et al1l200lL iBragg fc Kenvonl^OO^ . and the 

star NGC 884-2468 is a proposed candidate Be 
star ([Currie et al.l 120081 ). However, we do not see 
evidence of emission in our observations of these 
objects. As mentioned previously, the known Be 
star NGC 869-566 was initially observed by us in 
2005 and showed no sign of emission in its spec- 
trum. However, it has since developed increasingly 
stronger emission in our 2010 and 2012 observa- 
tions. For the other two transient Be stars (NGC 
869-49, NGC 884-1772), we do see emission in our 

spectra, however other authors have noted them 

i ; — 1 1 

in non-emission phases in the past (ISchildl 1196 
Slettebaklll985UBragg fc Kenvon|[2002l: iKeller et al 
20011 ). We are currently preparing a followup paper 
that discusses the Ha emission properties of all the 
Be stars in our sample. 

4. Comparison with Other Studies 



m 



Previous studies, su ch as lStrom et al] ([20051 ) and 
Huang fc Giesl ( 2006a ). have investigated some of 
our B star targets to determine the same basic stel- 
lar parameters we do, however their analyses em- 
ployed LTE atmospheric models or quantitative cor- 
rections to LTE model measurements to account for 
non-LTE effects. Between these two studies there 
are clear discrepancies in their measurements and 
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the conclusions each draws regarding the natal ro- 
tation rates and angular momentum evolution of the 
clusters' massive stars. An accurate determination 
of a given star's surface gravity is essential to the 
evaluation of stellar radius and evolutionary state. 
As we show here, non-LTE effects can contribute to 
significant errors in measurements of log g. 

Am ongst the 54 stars com mon to our sample and 
that of Huang fc Gied ( 2006a ) , we find some discrep- 
ancy between our respective results. As can be seen 
in Figure [6l there is very good agreement in our 
determinations of V sin i. This is to be expected 
since we both used LTE models to fit the He I and 
Mg II lines to measure V sin i. The differences at 
low V sin i may be due to a difference in spectral 
resolution between our respective datasets. 

However, differences in our measured tempera- 
tures for stars with T e g > 15000 K and in our mea- 
sured log g values are clearly apparent in Figures [7] 
and [51 respectively. Given that our general method- 
ology for measuring T R ff and log g is the same as 
Hua ng fe Gied (l2006ah . the source of these discrep- 
ancies lies part ially in our us e of dif ferent stellar 
models. While Huang fe Giesl ( 2006al) employ the 
LTE Kurucz ATLAS9 models (|Kurucd \l9M ), we 
use the more recently a vailable non-LTE TLUS TY 
BSTAR2006 models of lLanz fc Hubenvl (|2007l) for 
those stars in ou r sample with T P g > 15000 K. 



In their analysis iHuang fc Giesl (|2006al) acknowl- 
edge that their derived temperatures are likely to be 
slightly lower and gravities slightly higher compared 
to measurements derived from non-LTE model at- 
mospheres, as shown b y the comparative analysis of 
Lanz fc Hubeml (|2007l). 

Lanz fc Hubenvl ( 2007 ) find that in the non- 



LTE models the hydrogen Balmer lines tend to be 
broader and stronger due to the overpopulation of 
the n = 2 energy state, thus LTE models will yield 
overestimated surface gravities due to the altered 
the sh ape of the Balmer line wings. Przvbilla et al.l 
(|201lh compared LTE ATLAS9 models to non-LTE 
TLUSTY models for temperatures between 15000- 
35000 K. They found that non-LTE effects are sig- 
nificant above 22,000 K, affecting both the cores 
and wings of the Balmer lines. They find that LTE 
Balmer line profiles have equivalent widths up to 
30% lower than in non-LTE line profiles. This would 
cause a non-LTE model to find a higher temperature 
for the same observed line, or LTE models underes- 



find that LTE models of the H7 line may overesti- 
mate log g by up to 0.2 dex. 

While the expected temperature disagreement is 
opposite of the trend we find i n our c omparison of 
our work with Huang fc Gied ( 2006a ) results, sev- 
eral of our common sample stars with temperatures 
greater than 24,000K are Be stars, some of which 
we find to exhibit transient behavior. It is likely 
that emission has subtly filled in or otherwise altered 
the H7 line profile, which would res ult in the over- 
estimat ion of T e ff for these stars by IHuang fc Gied 
(2006a). We also find one proposed spectroscopic bi- 
nary among this common sample, so the H7 line pro- 
files may be further altered by variable line blending 
effects. For the remaining B-type stars in this re- 
gion it is possible that our temperature discrepancy 
is due to variable emission in unknown Be stars, un- 
resolved binaries, clumping in the hot stellar wind, 
or differences in the atomic species included in our 
respective atmosp heric models that affect h ydrogen 
stark broadening (jPrzvbilla fc Butlerl 120041 ) . 

With the quantified disagreement in T e $ and log 
g shown here, the anticipated effect of non-LTE 
atmospheres on the measurement r esults is more 
signific ant than initially assumed bv lHuang fc Gies 
(|2006al ). The effects of these overestimated tem- 
peratures and gravities also affect their determined 
stellar masses, their usage of log g as an indica- 
tor of stellar evolutionary status, and their deter- 
mined spin-down rates of B and Be stars. We find 
that clu ster members are mor e evolved than indi- 
cated bv lHuang fc Gied ( 2006a ). These B stars may 
spin down more slowl y tha n the rates observe d by 
Huang fc Giesl (|2006al) and lHuang eTail (|2010l ). 



We also com p are o ur results to those presented 
by IStrom et al. ( 20051 ) , who adopt the T e g values 
derive d from UBV photometry by ISlesnick et al 
(|2002n . In their study they find reasonable agree- 



ment b etween their measurem ents of V sin i and 
those of Huang fc Gied ( 2006a ). though their results 
are systematic a lly 5% smaller than the results of 
Huang fc Gied ( 2006a ). Comparing the results of 



the 26 stars common to our two samples, we find a 
similar agreement and systematic underestimation 
of Strom et al. 's V sin i values when compared to 
our measurements, as is expected giv en the excellent 
agreem ent of our results with those of lHuang fc Gies 
(|2006al ). U pon further compar ison of our results. 



timate the temperature. iPrzvbilla et al.l (|201l[ ) also 



we fin d that IStrom et alj (|2005l ) and ISlesnick et al 



(2002) have overestimated T e g for hotter stars as 
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well. We note that the two most discrepant stars 
are both Be stars, suggesting that their H/3 emission 
contaminates the B-band brightness used to derive 
T c ff. Thus M* and R± for the Be stars as shown in 
Figure |H] are also likely overestimated. 

Finally, we note that many of our Be stars 
we re found to be possib le sp ectroscopic b inarie s 
by iHuang fc Giesl (l2006ah an d IStrom et all f|2005h . 
Since IHuang fc Giesl (|2006al) did not present mea- 
surements of T e ff or log g for many of their spectro- 
scopic binaries, not all of our measurements could 
be directly compared. Their classification as bina- 
ries may be inaccurate due to variable emission in 
their spectral lines. Further monitoring of their ra- 
dial velocities as well as their emission will clarify 
their status. 

5. Conclusions and Future Work 

We have measured V sin i, T e g, log g po iar, M t , 
and i?* for 104 B-type and Be star members of 
NGC 869 and NGC 884 using spectroscopic mod- 
eling techniques and calculations from Stromgren 
photometry. Our determined values for V sin i 
are in good a g reemen t with the earlier results of 
IHuang fc Giesl ( 2006a ). though there is some dis- 
crepancy in our measured temperatures and sur- 
face gravities due to our use of the more recently 
available non-LTE BSTAR2006 stellar models of 



Lanz fc Hubenvl (|2007l ) and the possible contami- 
nation of Be stars and spectroscopic binaries. Bc- 
cause of the resultin g over-estimation of log g, 
Huang fc G ics (2006a) have underestimated the re- 



tention of initial angular momentum by the cluster 
members. 

We find that the cluster members are significantly 
more evolved than found by previous measurements. 
We also identify 8 transient Be stars in h and \ Per. 
The Be stars in these clusters are also rotating more 
slowly than expected based upon other young open 
clusters. Further monitoring of the massive stel- 
lar constituents of these clusters, and their rotation 
rates is well warranted. 

In a forthcoming paper we will examine our sam- 
ple of Be stars and their disk structures in greater 
detail. We will use the SEDs of B-type cluster mem- 
bers to compare our derived distances with previous 
measurements. Using our determined stellar param- 
eters, we will be able to separate and examine the 
stellar and disk contributions to the total Be star 



+ disk system flux through their observed SEDs. 
From these SEDs and multiple observations of Ha 
we will investigate the Be disk radii, masses, and 
longterm variability. 
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Fig. 1— Sample spectral line fits for NGC 869-90. 
Shown on the top is H7 and on the bottom is He I 
A4378. The solid line is our observed spectrum while 
the dashed line displays our model fit to the line, 
with the computed residual shown above, shifted 
for clarity. 
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Fig. 2. — Cumulative distribution function of V sin 
i for the Be stars (dashed line) and the normal B- 
type stars (dotted line) of both NGC 869 and NGC 
884. 
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Fig. 3. — T e ff measured for the B-type tempera- 
ture calibration stars from our work (diamonds) and 



from iNapiwotzki et al.l (1 1 993 | ) ( triangles ) compared 
to their calculated Tsaiona ( Balonal ll984). A linear 
fit of the two temperature scales (solid line) and the 
1:1 agreement (dotted line) are also shown. 
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Fig. 4. — Stromgren Co index and /3 magnitude are 
plotted with the cn — relations for luminosi ty class 
V and III stars ([Balona fc Shobbrooklll984t dotted 
lines). The B-type temperature calibration stars 
from this work (diamonds) are plotted with fifteen 
Be stars (filled diamonds) to demonstrate that the 
Be star (3 magnitudes are brightened due to the disk 
emission present. A best fit line for the calibration 
stars is also shown (solid line). 
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Fig. 5.— For both NGC 869 (top) and NGC 884 
(bottom), T e s and log g PO iar are plo t ted w ith the 
evolutionary tracks of ISchaller et al. | (|l992h - The 
zero age main sequence (ZAMS) mass of each evo- 
lutionary track is labeled along the bottom. Normal 
B-type stars are shown as open diamonds while Be 
stars are filled diamonds. 
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Fig. 6. — Compar ison of our V sin i m easure- 
ments with those of iHuang fc Giesl (|2006al ). NGC 
869 cluster members are shown as open diamonds, 
while NGC 884 members are shown as open trian- 
gles. Be stars are shown as filled diamonds and tri- 
angles. Spectroscopic binaries, as noted in Tables [2] 
and[3j are highlighted by double-sized symbols. 



Fig. 8. — Comparison of our log q measurements 
with those o f lHuang fc Giesl (|2006ah . in same format 
as Figure |U 
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Fig. 7. — Comparison of our T e s measurements 
with those o f lHuang fc Giel ()2006aj ) , in same format 
as Figure [6] 
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Fig. 9. — C omparison of our T e g measurements 
with t hose of Strom et al. ( 2005I) and lSlesnick et al 
(|2002j), in same format as Figure El 
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Table 2 

Physical Parameters of B-type Cluster Members 
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2119 


50 
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50 
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0.02 
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3 . 9 
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57 


3.84 


0.03 


4.07 
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0.1 
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7 
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0.10 


3.07 


21.0 


4.0 


22.2 
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5 
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0.02 
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8 


25500 
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3.37 


0.05 


3.40 


17.7 


0.8 


13.9 


1.9 


2392 


87 


17 


21200 


150 


3.75 


0.02 


3.81 


8.8 


0.1 


6.1 


1.3 


2520 


45 


10 


23200 


350 


3.67 


0.02 


3.68 


11.4 


0.3 


8.0 


1.7 


2540 


271 


23 


19800 


250 


3.53 


0.02 


3.82 


7.9 


0.2 


5.7 


1.1 


2555 


129 


6 


18250 


50 


4.05 


0.02 


4.13 


6.0 


0.1 


3.5 


1.3 


2605 


54 


27 


24600 
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3.58 


0.02 


3.60 


13.6 


0.5 


9.7 


1.9 


2622 


373 


35 


17200 


200 


3.38 


0.03 


3.91 


6.1 


0.1 


4.6 


0.8 


2716 a 


198 


7 


19650 


50 


3.75 


0.02 


3.92 


7.3 


0.1 


4.9 


1.2 


2729 


208 


17 


10410 


74 


3.83 


0.03 


4.01 


2.8 


0.1 


2.7 


0.1 


2907 


233 


52 


13700 


150 


3.86 


0.03 


4.07 


3.9 


0.1 


3.0 


0.4 



a Proposcd spectroscopic binary from Huans &l Gics (2006b) 

b Eclipsing binary fromJSouthworth ct al] l |2004|) 

c Proposed eclipsing binary from[Sacscn et alj t2OT0l) 



Table 3 

Physical Parameters of Be Star Cluster Members 



WEBDA V sin i AV sin i T cff AT cff log g A log g log g po i ar M, AM, fl, AR t Notes 
ID (km s" 1 ) (km s" 1 ) (K) (K) (dex) (dex) (dex) (M©) (M ) (Rq) (R e ) 
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49 
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21 


25400 


1790 


3.54 


0.23 


3.54 


15.4 


2.5 


11.0 
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No Emission Obscrved(l); 


























Emission Obscrvcd(2,7,8) 
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1115 


6 


25600 


600 


3.66 


0.05 


3.80 


12.7 


0.7 


7.4 


2.6 


No Emission Obscrvcd(2,8); 


























Emission Obscrvcd(4) 


309 


- 


- 


27700 


1790 


3.74 


0.23 


3.74 


15.8 


1.6 


8.8 


2.1 


Emission Obscrvcd(l,2,4,5,8) 


517 


178 


13 




















Emission Obscrvcd(4,6,7,8) 


566 a 


306 


25 


22100 


100 


3.25 


0.02 


3.63 


10.8 


0.1 


8.4 


1.2 


No Emission Obscrved(l,4,8); 


























Emission Obscrvcd(2,7,8) 
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22 


26700 


650 


3.91 


0.05 


3.96 


12.5 


0.8 


6 . 1 


3.2 


No Emission Obscrvcd(6,8); 


























Emission Obscrvcd(3) 
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205 


19 




















Emission Obscrvcd(4.6.7.8) 


847 


87 


11 


23000 


1790 


3.41 


0.24 


3.41 


13.7 


0.2 


12.1 


2.8 


Emission Obscrvod(2,4,6,7,8) 


992 






21700 


1790 


3.60 


0.24 


3.60 


10.7 


0.1 


8.6 


2.0 


Emission Obscrvcd(6,8) 


1057 






21300 


1790 


2.69 


0.25 


2.69 


29.5 


0.2 


40.8 


9.6 


Emission Obscrvcd(6.8) 


1161 






23400 


1790 


3.77 


0.24 


3.77 


11.1 


0.2 


7.2 


1.7 


Emission Obscrved(2,4,5,6,7,8) 


1261 


285 


78 


26400 


1790 


3.72 


0.23 


3.72 


14.3 


0.6 


8.6 


2.0 


Emission Observed(2,4,5,6,7,8) 


1268 


151 


12 


25400 


500 


3.46 


0.05 


3.59 


14.6 


0.6 


10.1 


2.2 


No Emission Obscrvcd(6,8); 


























Emission Obscrvcd(4) 


1278 


197 


12 


25100 


1790 


4.24 


0.24 


4.24 


9.9 


0.5 


3.9 


0.9 


Emission Obscrvcd(4.6,7,8) 


1282 






22700 


1790 


3.95 


0.24 


3.95 


9.2 


0.2 


5.3 


1.2 


Emission Observed(2,4,5,6,7,8) 


NGC 884 


























1702 a 


- 


- 


24400 


1790 


3.67 


0.24 


3.67 


12.6 


0.3 


8.6 


2.0 


Emission Observed(l,2,4,6,7,8) 


1772 


379 


28 




















No Emission Obscrvcd(4); 


























Emission Obscrvcd(6,8) 


1926 a 


106 


17 


27400 


1790 


3.89 


0.23 


3.89 


13.9 


1.2 


7.0 


1.6 


Emission Observed(2,4,6,7,8) 


1977 






20400 


1790 


4.27 


0.25 


4.27 


6.8 


0.1 


3.2 


0.7 


Emission Obscrvod(2,4,6,7,8) 


2088 a 






21900 


1790 


3.72 


0.24 


3.72 


10.0 


0.1 


7.2 


1.7 


Emission Obscrvod(2,4,5,6,7,8) 


2091 a 


236 


40 


18800 


1790 


3.97 


0.26 


3.97 


6.7 


0.1 


4.4 


1.0 


Emission Obscrvcd(4,6,7,8) 


2138 


153 


37 


24100 


1790 


3.53 


0.24 


3.53 


13.8 


0.2 


10.5 


2.5 


Emission Observed(l,2,4,5,7,8) 


2165 






24800 


1790 


3.79 


0.24 


3.79 


11.9 


0.3 


7.3 


1.7 


Emission Observed(2,4,5,6,7,8) 


2262 


265 


19 


24900 


1400 


3.51 


0.08 


3.75 


12.3 


1.4 


7.7 


2.2 


No Emission Obscrvcd(4,8); 


























Emission Obscrvcd(2,6) 


2284 






26500 


1790 


3.77 


0.23 


3.76 


14.0 


0.7 


8.1 


1.9 


Emission Observed(2,4,5,6,7,8) 


2468 


134 


46 


10070 


48 


3.88 


0.03 


4.00 


2.7 


0.1 


2.7 


0.1 


No Emission Obscrvcd(8); 


























Emission Obscrved(7) 


2563 


308 


74 


26100 


1790 


4.02 


0.23 


4.02 


11.6 


1.5 


5.5 


2.7 


Emission Observed(2,4,6,7,8) 


2771 






21800 


1790 


3.99 


0.24 


3.99 


8.5 


0.1 


4.9 


1.1 


Emission Obscrvcd(2,4,7,8) 



References: (1) Schild (1966), (2) Slcttcbak (1985), (3) Fabrcgat ct al. (1994), (4) Keller ct al. (2001), (5) Slcsnick ct al. (2002), (6) 
Bragg & Kcnyon (2002), (7) Be candidate with observed IR excesses in lCurrie et all i2008l) , (8) This work. 



^Proposed spectroscopic binary from Huang &: Gics (2006b) 



b Candidatc spectroscopic binary from Strom ct al. (2005) 



